Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region.

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

Contrasting patterns of nuclear and mtDNA diversity in Native American populations.

We report an integrated analysis of nuclear (autosomal, X- and Y-chromosome) short tandem repeat (STR) data and mtDNA D-loop sequences obtained in the same set of 22 Native populations from across the Americas. A north to south gradient of decreasing population diversity was observed, in agreement with a settlement of the Americas from the extreme northwest of the continent. This correlation is stronger with "least cost distances," which consider the coasts as facilitators of migration. Continent-wide estimates of population structure are highest for the Y-chromosome and lowest for the autosomes, consistent with the effective size of the different marker systems examined. Population differentiation is highest in East South America and lowest in Meso America and the Andean region. Regional analyses suggest a deviation from mutation-drift equilibrium consistent with population expansion in Meso America and the Andes and population contraction in Northwest and East South America. These data hint at an early divergence of Andean and non-Andean South Americans and at a contrasting demographic history for populations from these regions.

HLA-DR and HLA-DQ alleles in patients from the south of Brazil: markers for leprosy susceptibility and resistance.

BACKGROUND: Many epidemiological studies have shown that the genetic factors of the host play a role in the variability of clinical response to infection caused by M. leprae. With the purpose of identifying genes of susceptibility, the present study investigated the possible role of HLA-DRB1 and DQA1/DQB1 alleles in susceptibility to leprosy, and whether they account for the heterogeneity in immune responses observed following infection in a Southern Brazilian population. METHODS: One hundred and sixty-nine leprosy patients and 217 healthy controls were analyzed by polymerase chain reaction amplification and reverse hybridization with sequence-specific oligonucleotide probes and sequence-specific primers (One Lambda, CA, USA). RESULTS: There was a positive association of HLA-DRB1*16 (*1601 and *1602) with leprosy per se (7.3% vs. 3.2%, P = 0.01, OR = 2.52, CI = 1.26-5.01), in accord with previous serological studies, which showed DR2 as a marker of leprosy. Although, HLA-DQA1*05 frequency (29.8% vs. 20.9%, P = 0.0424, OR = 1.61, CI = 1.09-2.39) was higher in patients, and HLA-DQA1*02 (3.0% vs. 7.5%, P = 0.0392, OR = 0.39, CI = 0.16 - 0.95) and HLA-DQA1*04 (4.0% vs. 9.1%, P = 0.0314, OR = 0.42, CI = 0.19 - 0.93) frequencies lower, P-values were not significant after the Bonferroni's correction. Furthermore, HLA-DRB1*1601 (9.0% vs. 1.8%; P = 0.0016; OR = 5.81; CI = 2.05-16.46) was associated with susceptibility to borderline leprosy compared to control group, and while HLA-DRB1*08 (11.2% vs. 1.2%; P = 0.0037; OR = 12.00; CI = 1.51 - 95.12) was associated with susceptibility to lepromatous leprosy, when compared to tuberculoid leprosy, DRB1*04 was associated to protection. CONCLUSION: These data confirm the positive association of HLA-DR2 (DRB1*16) with leprosy per se, and the protector effect of DRB1*04 against lepromatous leprosy in Brazilian patients.

Otimização de metodologia PCR-SSP para identificação de polimorfismos genéticos de TNF e IL2 / Optimization of the PCR-SSP methodology in the identification of TNF and IL2 genetic polymorphisms

A análise de polimorfismos únicos de nucleotídeos (SNPs) de citocinas pode ser útil em estudos de frequências alélicas e genotípicas em populações saudáveis de diversas regiões, em estudos de associação com doenças infecciosas ou autoimunes, em estudos antropológicos e na evolução pós-transplante. Estes SNPs podem ser avaliados por diferentes métodos moleculares. O objetivo deste estudo foi aperfeiçoar uma metodologia PCR-SSP simples e rápida para a genotipagem de três SNPs de citocinas usando um único teste laboratorial. Para a identificação de IL2-330T/G e IL2+166G/T foram utilizados dois procedimentos na mesma genotipagem, cada um baseado no uso de quatro iniciadores. Para a detecção de TNF-238G/A foram utilizados dois iniciadores que amplificam a guanina e adenina na posição -238. Este estudo permitiu aperfeiçoar um método simples e rápido para identificar três SNPs de citocinas num único teste, podendo ser utilizado em qualquer laboratório de biologia molecular, como alternativa ao uso de kits de alto custo.

The analysis of cytokine single nucleotide polymorphisms (SNPs) can be useful in studies of allelic and genotypic frequencies in healthy populations from different regions of Brazil, in association studies of infectious or auto-immune diseases, in anthropological studies and in studies on post-transplant evolution. These SNPs can be assessed by different molecular methods. The objective of this study was to improve a simple and fast methodology, PCR-SSP, for the genotyping of three cytokine SNPs using a single laboratorial test. To identify IL2-330T/G and IL2+166G/T, two procedures were used in the same genotyping assay, each one based on the use of 4 primers. To detect TNF-238G/A, two primers were used that amplify guanine and adenine at position -238. This study enabled the improvement of a simple and fast method for identifying three cytokine SNPs in a single test, which can be adopted in any Molecular Biology laboratory as an alternative to the use of expensive kits.

Linkage disequilibrium patterns and genetic structure of Amerindian and non-Amerindian Brazilian populations revealed by long-range X-STR markers.

The extent of X-chromosome linkage disequilibrium (LD) was studied in a southern Brazilian population, and in a pool of samples from Amerindian populations. For this purpose, 11 microsatellites, located mostly in a Xq region comprising approximately 86 Mb was investigated. The lower Amerindian gene diversity associated with significant differences between the populations studied indicated population structure as the main cause for the higher LD values in the Amerindian pool. On the other hand, the LD levels of the non-Amerindian Brazilian sample, although less extensive than that of the Amerindians, were probably determined by admixture events. Our results indicated that different demographic histories have significant effects on LD levels of human populations, and provide a first approach to the X-chromosome ancestry of Amerindian and non-Amerindian Brazilian populations, being valuable for future studies involving mapping and population genetic studies.

Analysis of the CCR5 gene coding region diversity in five South American populations reveals two new non-synonymous alleles in Amerindians and high CCR5*D32 frequency in Euro-Brazilians.

The CC chemokine receptor 5 (CCR5) molecule is an important co-receptor for HIV. The effect of the CCR5*D32 allele in susceptibility to HIV infection and AIDS disease is well known. Other alleles than CCR5*D32 have not been analysed before, neither in Amerindians nor in the majority of the populations all over the world. We investigated the distribution of the CCR5 coding region alleles in South Brazil and noticed a high CCR5*D32 frequency in the Euro-Brazilian population of the Paraná State (9.3%), which is the highest thus far reported for Latin America. The D32 frequency is even higher among the Euro-Brazilian Mennonites (14.2%). This allele is uncommon in Afro-Brazilians (2.0%), rare in the Guarani Amerindians (0.4%) and absent in the Kaingang Amerindians and the Oriental-Brazilians. R223Q is common in the Oriental-Brazilians (7.7%) and R60S in the Afro-Brazilians (5.0%). A29S and L55Q present an impaired response to ß-chemokines and occurred in Afro- and Euro-Brazilians with cumulative frequencies of 4.4% and 2.7%, respectively. Two new non-synonymous alleles were found in Amerindians: C323F (g.3729G > T) in Guarani (1.4%) and Y68C (g.2964A > G) in Kaingang (10.3%). The functional characteristics of these alleles should be defined and considered in epidemiological investigations about HIV-1 infection and AIDS incidence in Amerindian populations.

Analysis of the CCR5 gene coding region diversity in five South American populations reveals two new non-synonymous alleles in Amerindians and high CCR5*D32 frequency in Euro-Brazilians

The CC chemokine receptor 5 (CCR5) molecule is an important co-receptor for HIV. The effect of the CCR5*D32 allele in susceptibility to HIV infection and AIDS disease is well known. Other alleles than CCR5*D32 have not been analysed before, neither in Amerindians nor in the majority of the populations all over the world. We investigated the distribution of the CCR5 coding region alleles in South Brazil and noticed a high CCR5*D32 frequency in the Euro-Brazilian population of the Paraná State (9.3 percent), which is the highest thus far reported for Latin America. The D32 frequency is even higher among the Euro-Brazilian Mennonites (14.2 percent). This allele is uncommon in Afro-Brazilians (2.0 percent), rare in the Guarani Amerindians (0.4 percent) and absent in the Kaingang Amerindians and the Oriental-Brazilians. R223Q is common in the Oriental-Brazilians (7.7 percent) and R60S in the Afro-Brazilians (5.0 percent). A29S and L55Q present an impaired response to beta-chemokines and occurred in Afro- and Euro-Brazilians with cumulative frequencies of 4.4 percent and 2.7 percent, respectively. Two new non-synonymous alleles were found in Amerindians: C323F (g.3729G > T) in Guarani (1.4 percent) and Y68C (g.2964A > G) in Kaingang (10.3 percent). The functional characteristics of these alleles should be defined and considered in epidemiological investigations about HIV-1 infection and AIDS incidence in Amerindian populations.

Two new mutations of the human BCHE gene (IVS3-14T>C and L574fsX576).

The genetic variation of human butyrylcholinesterase is associated with the majority of prolonged cases of apnea in patients submitted to the muscle relaxant succinylcholine. The present study reports two new mutations of the BCHE gene in 346 Euro-Brazilians: IVS3-14T>C found in five heterozygotes (allele frequency: 0.72+/-0.32%) and L574fsX576 found in one heterozygote (allele frequency: 0.14+/-0.14%). These two variants were not found in 85 Guarani Amerindians. It is not expected that the IVS3-14T>C mutation may interfere in the splicing process and that the mutation found in exon 4 (L574fsX576) may disturb BChE tetramerization and activity.

Geographic patterns of genome admixture in Latin American Mestizos.

The large and diverse population of Latin America is potentially a powerful resource for elucidating the genetic basis of complex traits through admixture mapping. However, no genome-wide characterization of admixture across Latin America has yet been attempted. Here, we report an analysis of admixture in thirteen Mestizo populations (i.e. in regions of mainly European and Native settlement) from seven countries in Latin America based on data for 678 autosomal and 29 X-chromosome microsatellites. We found extensive variation in Native American and European ancestry (and generally low levels of African ancestry) among populations and individuals, and evidence that admixture across Latin America has often involved predominantly European men and both Native and African women. An admixture analysis allowing for Native American population subdivision revealed a differentiation of the Native American ancestry amongst Mestizos. This observation is consistent with the genetic structure of pre-Columbian populations and with admixture having involved Natives from the area where the Mestizo examined are located. Our findings agree with available information on the demographic history of Latin America and have a number of implications for the design of association studies in population from the region.

Y-STR analysis in Brazilian and South Amerindian populations.

A sample of 203 Brazilian males from Rio Grande do Sul (RS), the Brazilian southernmost state, was typed for 11 Y-STR markers (DYS19, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS385, DYS437, DYS438, and DYS439). We also typed 42 individuals from two South Amerindian tribes (Kaingang and Guarani) to use the data as parental Amerindian contribution to our analyses. Gene and haplotypic diversities were estimated, with the South Amerindian samples showing smaller values for these parameters than Brazilians. To obtain a more comprehensive picture of the genetic structure of the Brazilian population as a whole, the Y-STR data from the RS sample was compared with those already published. No genetic substructuring was observed in the comparisons performed. Multidimensional scaling confirmed the proposed European source of most Y-chromosome Brazilian patrilineages.

Absence of the -116A variant of the butyrylcholinesterase BCHE gene in Guarani Amerindians from Mato Grosso do Sul

Butyrylcholinesterase (BChE; EC 3.1.1.8; Online Mendelian Inheritance in Man (OMIM) number 177400) is an enzyme found in many human tissues and encoded by the BCHE gene, of which 65 variants have been identified. In a recent study we found that the -116A variant of exon 1 of the BCHE gene was associated with lower mean BChE activity. The present study analyzed the -116 single nucleotide polymorphism (SNP) in 253 Guarani Amerindian Brazilians from the state of Mato Grosso do Sul (148 Guarani-Kaiowá, 83 Guarani-Ñandeva and 22 Kaiowá-Ñandeva descendants) and verified that they were all homozygotic for the -116G variant. A comparative analysis of the -116 site in nine vertebrate species indicated the -116A variant as the ancestral type. This is the first study of the -116 SNP in Amerindians and it is therefore difficult to infer whether or not the -116A variant was always absent from southern paleo-Amerindians or was present and then subsequently lost due to evolutionary factors.

Genetic variation and population structure in native Americans.

We examined genetic diversity and population structure in the American landmass using 678 autosomal microsatellite markers genotyped in 422 individuals representing 24 Native American populations sampled from North, Central, and South America. These data were analyzed jointly with similar data available in 54 other indigenous populations worldwide, including an additional five Native American groups. The Native American populations have lower genetic diversity and greater differentiation than populations from other continental regions. We observe gradients both of decreasing genetic diversity as a function of geographic distance from the Bering Strait and of decreasing genetic similarity to Siberians--signals of the southward dispersal of human populations from the northwestern tip of the Americas. We also observe evidence of: (1) a higher level of diversity and lower level of population structure in western South America compared to eastern South America, (2) a relative lack of differentiation between Mesoamerican and Andean populations, (3) a scenario in which coastal routes were easier for migrating peoples to traverse in comparison with inland routes, and (4) a partial agreement on a local scale between genetic similarity and the linguistic classification of populations. These findings offer new insights into the process of population dispersal and differentiation during the peopling of the Americas.

The beta-globin gene cluster distribution revisited-patterns in Native American populations.

New frequencies for the beta-globin gene cluster haplotypes are presented for the Aché (N = 82 individuals), Guarani (N = 76), and Kaingang (N = 54), three Native South American populations that live in an area between parallels 20 degrees S and 30 degrees S not covered by previous studies at this locus. The haplotype frequencies obtained for the three populations are within the interval observed for 28 other Native American populations. The Aché show much less haplotypes (five) than the other two populations (9-10), the haplotype prevalences being more similar to those of the Guarani than to the Kaingang. The Native American total heterozygosity was about half (0.41) that obtained for the African populations (0.71), but was not much different from those obtained for other continents. A geographical pattern was disclosed in South America by mapping the frequencies of the most common haplotype (haplotype 2), and by means of spatial correlation analysis. The analysis of molecular variance (AMOVA) and pairwise F(ST) data suggest three distinct sectors for the genetic landscape of Native South America: the Andes, the Center/Southeast region, and the Amazon.

Butyrylcholinesterase genetic variability in Guarani Amerindians from the Brazilian state of Mato Grosso do Sul

Human butyrylcholinesterase (BChE; EC 3.1.1.8) is a polymorphic enzyme coded by the BCHE gene (3q26.1-q26.2) while the CHE2 gene (2q33-q35) determines a still not characterized substance that forms a complex with BChE (C5), being the CHE2 C5+ and CHE2 C5- phenotypes detected in electrophoresis. The present study investigated BCHE and CHE2 variability and the BChE activity of Brazilian Guarani Amerindians from the Kaiowá and Ñandeva sub-groups living in several indigenous territories in the Brazilian state of Mato Grosso do Sul. The frequency of the BCHE exon 2 D70G (A) allele was 0.60 percent ± 0.35 percent while that of the BCHE exon 2 G390V (F-2) allele, never before screened in Amerindians, was 8.82 percent ± 1.35 percent. This is the first time that the BCHE gene exon 4 A539T (K) allele has been surveyed in Brazilian Amerindians where it was found at a frequency of 3.69 percent ± 0.85 percent, similar to that found in Chilean Mapuche Amerindians. The BCHE gene variability seen in this survey differs from that of non-isolated populations in respect to both A539T and G390V allele frequency. The CHE2 C5+ phenotype frequency was 14.40 percent ± 2.22 percent and falls within the range of that found for other Brazilian Amerindian samples.

Human T - cell lymphotropic virus type II in Guaraní Indians, Southern Brazil

O vírus linfotrópico de células T-humanas do tipo II (HTLV-II) é identificado em muitos grupos de ameríndios. No Brasil, tem sido encontrado em indivíduos da população urbana, bem como em índios oriundos da região Amazônica. Os Índios Guaraní, do Sul do país, foram investigados para infecção por HTLV-I/II. Três indivíduos, oriundos de uma amostra de 52 índios, demonstraram sororeatividade para HTLV-II (ensaio imunoenzimático e Western blot). Este estudo preliminar foi o primeiro a identificar a presença de infecção por HTLV-II em ameríndios do Sul do Brasil.

Human T-cell lymphotropic virus type II in Guaraní Indians, Southern Brazil.

Human T-cell lymphotropic virus type II (HTLV-II) is found in many New World Indian groups on the American continent. In Brazil, HTLV-II has been found among urban residents and Indians in the Amazon region, in the North. Guaraní Indians in the South of Brazil were studied for HTLV-I/II infection. Among 52 individuals, three (5.76%) showed positive anti-HTLV-II antibodies (enzyme-linked immunosorbent assay and Western blot). This preliminary report is the first seroepidemiological study showing HTLV-II infection among Indians in the South of Brazil.

Geography influences microsatellite polymorphism diversity in Amerindians.

Data related to 15 short tandem repeat polymorphisms (STRPs) are reported for four South American Indian populations, and integrated with previous Brazilian Indian results. Overall heterozygosities varied significantly among groups (Kruskal-Wallis test, P = 0.002). The lowest levels of heterozygosity were observed in the Ache, Ayoreo, and Surui, an expected finding considering their isolation and ethnohistory. Genetic distance and gene diversity analyses suggested that geography was a good predictor of genetic affinity among these Native Americans. New evidence from this study supports the hypothesis that the Ache population descends from a Ge group that preceded the Guarani colonization of Paraguay.

Y-chromosome evidence for differing ancient demographic histories in the Americas.

To scrutinize the male ancestry of extant Native American populations, we examined eight biallelic and six microsatellite polymorphisms from the nonrecombining portion of the Y chromosome, in 438 individuals from 24 Native American populations (1 Na Dené and 23 South Amerinds) and in 404 Mongolians. One of the biallelic markers typed is a recently identified mutation (M242) characterizing a novel founder Native American haplogroup. The distribution, relatedness, and diversity of Y lineages in Native Americans indicate a differentiated male ancestry for populations from North and South America, strongly supporting a diverse demographic history for populations from these areas. These data are consistent with the occurrence of two major male migrations from southern/central Siberia to the Americas (with the second migration being restricted to North America) and a shared ancestry in central Asia for some of the initial migrants to Europe and the Americas. The microsatellite diversity and distribution of a Y lineage specific to South America (Q-M19) indicates that certain Amerind populations have been isolated since the initial colonization of the region, suggesting an early onset for tribalization of Native Americans. Age estimates based on Y-chromosome microsatellite diversity place the initial settlement of the American continent at approximately 14,000 years ago, in relative agreement with the age of well-established archaeological evidence.

Polymorphisms of CYP1a1, CYP2e1, GSTM1, GSTT1, and TP53 genes in Amerindians.

Polymorphisms at the TP53, cytochrome P-450 (CYP), and glutathione S-transferase (GST) genes are related to cancer susceptibility and present high diversity in allele frequencies among ethnic groups. This study concerns the CYP2E1, GSTM1, and GSTT1 polymorphisms in seven Amerindian populations (Xavante, Guarani, Aché, Wai Wai, Zoró, Surui, and Gavião). Polymorphic sites at CYP1A1 and TP53 were also studied in the Aché and Guarani tribes and compared with previous results about these systems already obtained in the other populations. The CYP2E1*5B haplotype showed, respectively, the highest and the lowest frequencies already observed in human groups. High frequencies of CYP1A1*2A and CYP1A1*2C alleles and mostly low values of GSTM1*0/*0 and GSTT1*0/*0 genotypes were observed. These data may be interpreted as being due to genetic drift or selection for these high-frequency CYP1A1 alleles and against GST null genotypes during America's colonization. Intrapopulation diversity varied from 0.19 (Guarani) to 0.38 (Surui), and 90% of the total diversity was due to the variability within populations. The relationships between these Amerindians and with other ethnic groups were evaluated based on D(A) distances and the neighbor-joining method. Low correlation was observed between genetic relationships and geographic distances or linguistic groups. In the TP53 comparison with other ethnic groups, Amerindians clustered together and then joined Chinese populations. The cluster analysis seems to indicate that the Aché tribe might descend from a Gê group that could have first colonized that Paraguayan region, but had also assimilated some amount of the Guarani gene pool, maybe through intertribal admixture.